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Abstract 

We study dynamical aspects of three-dimensional gonihedric spins by using Monte-Carlo 
methods. The interest of this family of models (parametrized by one self-avoidance pa- 
rameter k) lies in their capability to show remarkably slow dynamics and seemingly glassy 
behaviour below a certain temperature T g without the need of introducing disorder of any 
kind. We consider first a hamiltonian that takes into account only a four-spin term (k = 0), 
where a first order phase transition is well established. By studying the relaxation prop- 
erties at low temperatures we confirm that the model exhibits two distinct regimes. For 
T g < T < T c , with long lived metastability and a supercooled phase, the approach to equi- 
librium is well described by a stretched exponential. For T < T g the dynamics appears 
to be logarithmic. We provide an accurate determination of T g . We also determine the 
evolution of particularly long lived configurations. Next, we consider the case k = 1, where 
the plaquette term is absent and the gonihedric action consists in a ferromagnetic Ising 
with fine-tuned next-to- nearest neighbour interactions. This model exhibits a second order 
phase transition. The consideration of the relaxation time for configurations in the cold 
phase reveals the presence of slow dynamics and glassy behaviour for any T < T c . Type II 
aging features are exhibited by this model. 
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1 Introduction 



Glassy systems are very common in nature, yet not quite well understood. Lattice models 
may serve as good candidates to describe some properties of these systems. In recent years 
some interest has been raised by remarkably simple Ising spin systems that originated from 
the consideration of a model of random surfaces in the context of string theory p[. The 
version of this model in a discretized space defines the so-called gonihedric spin model which 
consists, in general, of an Ising model with finely tuned nearest, next-to-nearest neighbour and 
plaquette interactions. The relation among the couplings of the hamiltonian depends also on the 
dimensionality of the system. The geometric origins of the model show up in a remarkable simple 
way of writing the energy of a given configuration: the surfaces corresponding to the interfaces 
between up and down spins are weighed by E = n 2 + Kn 4 , where n 2 is the number of edges of such 
an interface and is the number of four plaquettes that share a common link. The parameter 
k can thus be interpreted as an indicator of the self-avoidance of the model. Notice that there is 
no microscopic surface tension. 

Up to now there has been a considerable amount of numerical work on the three-dimensional 
case, to which we shall refer in the following and also some peliminiary results in four-dimensional 
case ||. In two dimensions, the model with k — is actually trivial (no phase transition |4j]) but 
the solution with k ^ is unknown. 

For the three dimensional case which is of our concern in this paper the Hamiltonian of the 
model takes the form 

H(a) = -2/t^ a r -CT r - +cJ +7;J2 ^r+a+Z? g"" °r°?+aOf +Si+ 0O?+p, (1) 

where a and (3 are lattice unit vectors. The model is defined on an cubic lattice. The system 
exhibits a very high degree of symmetry due to the particular ratio of the couplings. This 
symmetry implies that there is no cost by the flipping of any plane of spins. This results in a 
highly degenerate ground state. Note in passing that this last feature is common in every glassy 
system. 

Some numerical evidence regarding dynamical and equilibrium properties of the hamiltonian 
([!]) has been accumulated in the past, particularly for n = 0, but also for k ^ 0. These studies 
revealed many interesting features and provided evidence that a glassy phase is present in the 
phase diagram 0. First of all, it is known that the model has a first order phase transition at a 
temperature T c where the solid to liquid transition is present I|. There is also good evidence that 
a dynamical transition exists at T g < T c which seems to mark the onset of the glassy behaviour 
0. The existence of T g and the study of the dynamical properties of the system above and below 
that value have been considered in f)|-||13|. Preliminary results for k ^ were given in |14| in 



what respects the equilibrium properties, and in |15j in what respects the dynamical properties 
of the system . 

In section 2 we provide a qualitative view of the slow dynamics behaviour of the cold phase 
by considering the relaxation properties of the model and we estimate with good accuracy T g 
by measuring the spin-spin autocorrelation functions. We show that there is a dramatic change 
in the behaviour of this function above and below T g . In section 3 we study a different case 
of the spin gonihedric action by taking k = 1 in ([!]), apparently a much simpler system (the 

^For the interesting case of the four-spin model where randomly distributed couplings are considered, see (pf). 
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plaquette term in the hamiltonian is absent for this value). For that value of the parameter k 
the system is described by nearest and next-to-nearest interactions. We confirm, by studying 
the energy susceptibility, that the system has a second order phase transition f\ . Furthermore 
we find slow dynamics behaviour anywhere in the cold phase (below the critical temperature). 
By studying spin-spin autocorrelation function and the overlap spin distribution function, we 



provide evidence that the system exhibits type II aging jl7| , which is a feature of the glassy 
systems. 

The physical interest of this model is twofold. On the formal side, the model is of such 
simplicity that a theoretical understanding of the mechanisms underlying slow dynamics and 
glassy behaviour appears possible. On a more practical side, it would be extremely interesting to 
be able to understand and produce magnetic materials and coatings with such finely tuned (or 
approximately so) couplings. The extremely long relaxation times would make them very robust 
agains thermal noise and fluctuations, yet encoding information there would be as simple as in a 



normal magnetic material. This possibility has been suggested in fL8| . While the plaquette term 
seems hard to imitate in real materials, the fact that many of the interesting features persist for 
k = 1 makes perhaps such possibility less remote. 



2 Four— spin interaction (k = 0) 

In this section we will study the case k = 0, i.e, a spin model with only four-spin (plaquette) 
interaction. In this case the hamiltonian (|1|) takes the form: 

H = \ <7^+S(T? +d £ + 0V? + (2) 

r,S,/3 

This form of the interaction leads to a highly degenerated ground state. Flipping every spin 
in any plane of the cubic lattice implies invariant ground state energy. Taking this symmetry 
into account, the degenerancy of the ground state is equal to 2 3L due to the 3L diferent planes 
in a cubic lattice. This degeneracy survives even at T ^ 0. 

It is well known that this model exhibits a first order phase transition at T c = 1.95 along 
with a dynamical transition at T g ~ 1.7, which is the temperature where the glassy behaviour 
shows up [f|, [|HJ, |pT| . Even though our main interest is the study of the glassy characteristics 
by looking at the relaxation as well as the autocorrelation of the order parameters (to be defined 
below) in the glassy phase, the region T g < T < T c is interesting as well. In this region, 
numerical simulations clearly indicate the presence of metastability. This is exemplified by the 
result presented in Fig. |l]a where the time evolution of an 8 3 volume at a temperature value 
just below T c is shown This figure also gives clear evidence that there is a first order phase 
transition. Fig. [l]b consists of two different curves. One of them corresponds to a heating process 
starting from an initial ordered congiguration. The other one describes the result of quenching 
a random initial configuration for each one of the temperature values shown. In both cases the 
simulation has been carried out using the Metropolis algorithm on a 24 3 volume by performing 10 4 
measurements at each temperature value. As we will see in a while, the approach to equilibrium 

"Note that from some preliminiary results referred in (Q) the first order transition present for n = gets 
weaker and possibly becomes second order at n ~ 0.5. 

**This simulation refers to the system expressed in dual variables [^| and it has been carried out using a cluster 
algorithm. 
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V«24 A 3 




in the region T g < T < T c is non-standard and it is well described by a stretched exponential, 
instead of a simple exponential. For T < T g one immediately sees that the results from the 
quench of the random configuration after 10 4 MC steps differ from those obtained starting from 
an ordered configuration (any of the 2 3L vacua) 0. The difference appears to be constant all 
over this region (for a fixed number of thermalization steps). This clearly hints to the coexistence 
of two different dynamics. Initially fast dynamics quickly brings an initial configuration which is 
badly out of equilibrium to some sort of approximate equilibrium. At that point slow dynamics 
takes over and the evolution of the system is considerably freezed. 

Before getting into the more quantitative aspects of these results, it is perhaps interesting to 
turn to one of our motivations, namely to test whether the appearence of slow dynamics makes 
the transition between two approximate ground states so slow as to make a given configura- 
tion virtually indestructible by thermal fluctuations, thus providing a convenient way of storing 
information. 



To this end, we simulate the system on a cubic lattice and we use a Glauber algorithm |JT9 
which is assumed to provide a good approximation to the thermal mechanism of fluctuations. 
Metropolis or heat bath give very similar results. We look at the decay of an artificial initial 
configuration consisted of an inner volume with a chessboard-like arrangement of the spins, 
which is one of the ground states, while the outer volume spins are fixed at +1 i.e, they form 
a ferromagnetic ground state. We use fixed boundary conditions in order to make the system 
decay to a ferromagnetic ground state. Our measured quantity is the number of the minus spins, 
N-(t), in terms of the computing time divided by the number of spins at time t = 0, AL(0). The 
quantity N-(t) is clearly related with the magnetization, M, since: M « N+ + N- = iV t0 i — 2iV_. 
Two examples of our results are given in Fig. 0a and Fig. 0b which have been produced from 
simulations over two lattice volumes namely 16 3 and 30 3 which at time t — have enclosed 8 3 
and 20 3 lattice volumes respectively with a chessboard-like arrangement for the spins. 

The two curves in each sub figure have been generated after averaging over a sample of 50-100 
copies where each one evolves starting from the same initial configuration, in order to reduce 
the noise. Two cases are shown. One is for temperature T = 1.83, in the supercooled phase, 
while the other one, T = 1.0, lies in the glassy phase . The difference on the relaxation time 
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(a) 



(b) 



Figure 2: The relaxation for two diferent temperatures, T = 1.0 and T = 1.83: (a) 16 3 lattice 
volume with an initial chessboard-like configuration in a 8 3 volume, (b) 30 3 lattice volume with 
an initial chessboard-like configuration in a 20 3 volume. 

between the two temperatures is obvious. For the T = 1.83 case the system reaches a stable 
value relatively fast for both lattice volumes, but this is not the case for T = 1.0. Note that the 
slope of the curve keeps taking non-zero value even for remarkably large times. This is perhaps 
clearer in Fig. |3| where long runs are depicted for T = 1.0, confirming that we are in presence of 
slow dynamics. This behaviour can be seen in a more apparent way in Fig. [| where the results 
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Figure 3: The long time evolution at T=1.0 for inner chessboard-like configuration in a 8 3 (a) 
and 20 3 volume (b). 

of Fig. [| are presented in a logarithmic time scale. In Fig. [| the logarithmic decay is present 
clearly enough. 

The very distinct dynamics between the glassy phase and the evolution of the supercooled 
phase in the metastability region T g < T < T c can be seen in a more definite way in Fig. |5]. 
In this plot the results for four temperatures are presented for the bigger volume used i.e. 30 3 . 



4 



100 



1000 



10000 



100000 



t (MC) 



Figure 4: The Fig. |3] results as plotted in a logarithmic time scale. (8 3 and 20 3 initial chessboard- 
like configurations for total volumes 16 3 and 30 3 respectively.) 



All of them correspond to a random starting configuration. For T = 1.83, which lies in the 
metastability region the evolution seems very fast all the way to the equilibrium value. On the 
contrary for the other three temperatures, a fast evolution is initially observed followed by a very 
slow one which persists up to very long times. 

These results are very suggestive and indeed show that the system finds very difficult to 
overcome dynamical energy barriers that are created along the evolution and this is undoubtedly 
the reason for the slow dynamics. Recall that the energy of the model is concentrated on the edges; 
the system has vanishing microscopic surface tension. To reduce the volume of the excitation 
with local moves, the total edge length must temporarily increase by a substantial amount. That 
makes excitations such as the one we have been analysing virtually stable. 

In spite of the unambiguity of the previous results, in order to determine some properties 
associated to the very slow dynamics observed at low temperatures we shall proceed to analysing 
several dynamical correlators. 

We shall first study the spin-spin autocorrelation function for temperatures lying in the 
metastability region after a random start, thus forcing the system to be in a supercooled phase. 
Its definition is given by 



The brackets mean that we take the average value over copies starting from a random config- 
uration (200-400 in our case). The waiting time t w is the time for which the system is being 
thermalized before taking the measurements at subsequent times denoted by t. The waiting time 
t w is taken to be about 300 by noticing that from that value on the resultant values for C sp i n are 
quite identical as long as the temperature lies in the metastability region. In Fig. ^| we present 
some of our results for C spin in a 40 3 lattice volume and for four temperature values. The fittings 
are stretched exponentials of the form ae~^/ T ) . We denote by r the relaxation time. In all 
cases we found 0.60 < b < 0.80 with an error smaller than 0.004. Since the fits seem quite good 
we go on and make a plot of the resulting values of r versus the temperature T. The result is 




(3) 
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Figure 5: For T = 1.83 the system is in the supercooled region and exhibits fast relaxation. For 
the three other temperatures values which stand in the cold phase the slow decay is obvious when 
it is presented in logarithmic time scale. 



shown in Fig. [7] where the corresponding fit to a function having the form, const/(T — T g ) c , is 
quite good and leads to the prediction T g = 1.698(1) with c = 0.41. At T g the autocorrelation 
time is expected to diverge because of the onset of the slow dynamics which turns the stretched 
exponential behaviour into a power law (with a small exponent) or a logarithm. 




Figure 6: Spin-spin autocorrelation function versus time. The fittings to the points are of 
stretched exponential type. 



Although that result is in good agreement with previous simulations |§, |TTJ the method 
relying on the stretched exponential fit may prove to be too risky for an exact prediction of 
T g due to ambiguities in the fitting process and perhaps is not trustable to that accuracy. As 
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Figure 7: r versus T. The fitting procedure leads to divergent r for T g = 1 .698. 



an alternative and a cross check, for two different temperature values, namely T = f .720 and 
T = 1.695 we show the behaviour of the C sp i n for two very different values of the waiting time, 
t w = 300, 2000. One sees that for the higher temperature the behaviour of the C sp i n is identical 
but for the smaller one there is a strong dependence on t w showing that for larger waiting times 
the system exhibits different, and in particular, slow dynamics behaviour which is a strong signal 
that the system has passed to the glassy phase. Hence we are able to estimate the value of T g 
to lie in between the interval 1.695 < T g < 1.720 which indeed agrees also with the previously 
obtained value. 



t(MC) 



T=1.720, tw=300 
T=1.720, tw=2000 
T=1.695, tw=300 
T=1.695, tw=2000 




Figure 8: C sp i n versus time for two values for the waiting time t w . For the system being in the 
glassy phase we see the dependence on t w . 
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3 Next and next-to-nearest neighbour interaction (k = 1) 



As we have indicated the gonihedric spin model is actually a family of models which are parametrized 
by a real quantity k. All the members of this family share the common feature of having their 
spin interfaces weighed with the total edge length. The value of k simply indicates the degree of 
self-avoidance of the surfaces. The plaquette model whose dynamical properties we just discussed 
corresponds to k — 0. An interesting member of this family is given by k — 1. The degeneracy 
of the ground state in this case is 3 • 2 L . 
In this case the Hamiltonian is given by 

H = -2 a r a r+a + \ ( 4 ) 

As we see, the spin plaquette term has disappeared. From a practical point of view this model 
may be particularly interesting as the plaquette term is obviously hard to get in real materials. 
From the standpoint of a spin system this is just a model with nearest and next-to-nearest 
neighbour interactions (though a finely tuned one). 

We shall proceed to studying the dynamical properties of the system by using Monte-Carlo 
methods. In performing the simulations we used the Metropolis algorithm for several lattice 
volumes, namely 10 3 , 16 3 , 20 3 , 24 3 , 30 3 , 40 3 and 46 3 . We imposed either fixed or periodic 
boundary conditions depending on the kind of the measurement that was performed. During the 
presentation of our results it will be explicitly mentioned when we used one or the other. 

We begin by giving the behaviour for the susceptibility of the energy with the lattice volume 
which is given by 

S(E) = V((E 2 ) - (E 2 )) 

We denote by V the lattice volume and by the symbol (), the average value over sweeps. In Fig. 
the susceptibility for the system energy as a function of the temperature for four lattice volumes 
is depicted. For every point in the figure we performed 10 5 thermalization sweeps followed by 
more than 10 5 measurements, using periodic boundary conditions and starting from the ordered 
configuration. The peak for each volume clearly increases with it although with an exponent 
less than one, which is a signal for a second order phase transition. Also, the positions of the 
peaks for the bigger volumes are seen to concetrate around the value T c = 2.329 which is the 
pseudocritical value for the up to the volumes we used. It should be mentioned that while the 
second order character of the phase transition was already mentioned in however our 

prediction for T c gives a somewhat smaller value. 

Having estimated the critical temperature value for the system we want now to study the 
dynamic behaviour at low temperatures, i.e. below T c . To this end we use first the same method 
as for the k = case, described in the previous section, which consists of choosing an initial 
chessboard-like configuration defined in a cubic subvolume of the system. We impose fixed 
boundary conditions and we study the behaviour of the dynamical quantity 



E(t) - E eq 
E(0) - E, 



(5) 



eq 



with time for various temperatures in the cold phase and for different sizes of the inner chessboard- 
like volume. In Eq.(|5|) E(t) and E(0) denote the energy of the system at time t and just after the 
first sweep performed, respectively and E eq is the system energy when an ordered configuration 
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Figure 9: S(E) versus T behaviour with the lattice volume. 



is taken as the initial one. In Fig. |T(1| a we give an example for lattice volume 16 3 . We fix the 
temperature at the value T = 1.1 and we are interested in the behaviour as the dimension of 
the inner chessboard-like volume increases from L = 2 to L = 8. Each curve has been produced 
after averaging several hundreds of repetitions starting from the same initial configuration in 
order for the noise to be reduced. The results show a dramatic increase of the relaxation time 
with increasing system size. For the maximum value of L used, the observable (0) shows very 
slow variation with the time and it follows a power law behaviour with estimated exponent equal 
to 0.18, so a logarithmic behaviour cannot be really excluded. 

We can reach similar conclusions by studying the above relaxation quantity for decreasing 
values of temperature at fixed L — 3. In Fig. |T0"| b the corresponding results are depicted for a 
16 3 total volume. In particular for the value T = 0.2 the dynamics is so slow that it seems really 
hard for the system to get to the ground state for accessible computing times. 

Next we shall move to a more quantitative study following the same lines of the previous 
section. We look for aging features by considering the autocorrelation spin-spin function defined 
by (|3]). Brackets indicate that we take the average value over 100-200 trials starting from a 
random configuration. t w is the time for which the system is being thermalized before taking 
the measurements at subsequent times denoted by t. We follow the same lines of analysis as in 



TT| which concerns the k = case. In Fig. |TT] we show an example of the behaviour of the 
autocorrelation function (^) for the temperature value T=2.2 and for a 24 3 lattice volume. Six 
different curves for the C sp i n are depicted each corresponding to a different value of t w . The very 
slow evolution of the autocorrelation function combined with the clear dependence on the time 
t w give strong evidence in favour of aging. In order to determine the aging (whether I or II) we 
consider the overlap function |T7|| : 

Q(t + t w ,t + t w ) = ^(Y,^\t + t w )ap\t + t w )) (6) 

r 

The measuring process for the above consists of the relaxation of the system for time t w 
starting from the random configuration. At that moment we take two copies and each 
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Figure 10: The relaxation behaviour of the initial configuration as seen by the study of quantity 
defined in ([5]): a) for various inner volume sizes at fixed temperature (from top to bottom, 
L = 2,3,4,6,8) and b) for fixed inner volume size with decreasing temperature (from top to 
bottom T = 0.2,0.5,0.9,1.1). 
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Figure 11: The spin-spin autocorrelation function as depicted for different waiting times t w (from 
top to bottom t w = 1000, 500, 200, 100, 50, 10). 




evolving independently for time t. In Fig. [12] the behaviour of Q(t + t w ,t + t w ) in terms of 
C sp i n (t, t w ) for the same values for the volume and temperatute as in Fig. |Tl| and for four values 
of the waiting time t w is shown. It clearly follows that as C spin (t, t w ) — > 0, Q(t + t w ,t + t w ) — > 0. 
The fact that the two system copies are moving independently since they show zero overlapping 
indicates strong evidence for type II aging. 
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Figure 12: The Q(t + t w ,t + t w ) — > behaviour as C spin {t,t w ) — > is clearly seen for the four 
values of t w . 



4 Conclusions 

At this point we would like to summarize our main conclusions. 

We have carried out a rather complete analysis of the dynamical properties of the gonihedric 
spin model. The model naturally exhibits large potential barriers that are dynamically generated 
along the evolution process when trying to adjust to the environment. This is due to the fact 
that the area term of the interfaces is completely subleading and the dynamics is driven by the 
total edge length of the excitations. This creates basins of stability from where it is virtually 
impossible to kick out the system by local thermal fluctuations. 

Several algorithms have been used: Glauber dynamics, heat bath and standard Metropolis. 
Modulo an overall rescaling of Monte-Carlo time, the results are fully equivalent. 

In this paper we have studied numerically two cases of the gonihedric spin model. For the 
k = case we give evidence of a glassy transition and we estimate the temperature value T g at 
which the glassy phase arises being in the interval [1.695,1.720]. We show that the dynamics 
is most likely logarithmic for T < T g . In the mestatability region, above T g but below the true 
thermodynamical transition T c (where the system has a first order transition) the approach to 
equilibrium is well described by a stretched exponential, whose characteristic time scale diverges 
as T -> T+. 

We confirm that for the case k — 1 the system exhibits a second order phase transition. The 
critical properties of this transition have not been elucidated yet; a tentative determination of the 



exponents using finite size scaling carried out in |14| gives a rather non-standard set of values. 
Work on this is in progress, but the problem is agravated by the presence of slow dynamics. 
Indeed we have detected in the cold phase of this system very slow dynamical behaviour, quite 
similar to what takes place with k = for T < T g . This dynamical behaviour has all the features 
of type II aging and is again compatible with logarithmic evolution. 

Possible applications of these simple, but interesting, spin models are discussed. 



11 



Acknowledgements 



D.E. thanks D. Johnston and A. Lipowski for discussions and the hospitality of the Department of 
Mathematics of Heriot-Watt University. P.D. is grateful to G. Koutsoumbas for useful discussions 
and acknowledges support from "EUROGRID- Discrete random geometries: from solid state 
physics to quantum gravity" (HPRN-CT-1999-00161). The support from grant MCyT FPA 
2001-3598 and CIRIT 200156R-00065 is also acknowledged. Some of this work was carried out 
at CESCA. 



References 

[1] R.V: Ambartzumian, G.S. Sukasian, G.K. Savvidy and K.G. Savvidy, Phys. Lett. B275 
(1992) 99; G.K. Savvidy and K.G. Savvidy Int. J. Mod. Phys. A8 (1993) 3393; G.K. 
Savvidy and K.G. Savvidy, Mod. Phys. Lett. A8 (1993) 2963. 

G.K. Savvidy and F.J. Wegner, Nucl. Phys. B413 (1994) 605; G.K. Savvidy and K.G. 
Savvidy, Phys. Lett. B324 (1994) 72; G.K. Savvidy, K.G. Savvidy and F.J. Wegner, Nucl. 
Phys. B443 (1995) 565; J. Ambjorn, G. Koutsoumbas and G.K. Savvidy, Europhys.Lett. 
46 (1999) 319 ( |cond-mat/9810271|) . 



G. Koutsoumbas, G.K. Savvidy and K.G. Savvidy, Phys.Lett. B410 (1997) 241 ( [hep - - 
th/9706173|) . 



G.K. Bathas, E. Floratos, G.K. Savvidy and K.G. Savvidy, Modern Physics Letters A10 
(1995) 2695 (|hep-th/9504054j ). 

D. Alvarez, S. Franz, F. Ritort, Phys. Rev. B54, 9756 (1996) (|cond-mat/9603200| ). 

D.Espriu, M. Baig, D.A.Johnston, R.K.P.C.Malmini, J.Phys. A30 (1997) 405 ([Eepj 
lat/9607002|) . 

G. Koutsoumbas and G.K. Savvidy, Three-dimensional gonihedric spin system ( |cond- 
mat/0111590D . 



E.N.M. Cirillo, G. Gonella, DA. Johnston and A. Pelizzola, Phys. Lett. A226 (1997) 59. 
A. Lipowski, J. Phys. A: Math. Gen. 30 (1997) 7365. 

A. Lipowski and D. Johnston, Glassy transition and metastability in four-spin Ising model 
(|cond-mat/9812098|) . 



M.R. Swift, H. Bokil, R.D.M. Travasso and A.J. Bray, Phys.Rev. B62 11494 (2000) flcond- 
mat/0003384T ). 



A. Lipowski and D. Johnston, Phys.Rev. E61, 6375 (2000) flcon3"-mat/9910370D . 

A. Lipowski and D. Johnston, Phys.Rev. E64, 041605 (2001)( |c^Hcrmat/0l05602l ). 

M. Baig, D.Espriu, D. Johnston and R.P.K.C. Malmini, J.Phys. A30 (1997) 7695 Qnep- 
lat/9703008|) . 



12 



[15] A.Lipowski, D.Johnston and D.Espriu, Phys.Rev. E62 (2000) 3404 (|cond-mat/000446q ). 
[16] D.A. Johnston and R.P.K. Malmini, Phys.Lett. B378 (1996) 87, Qhep-lat/9508026; ). 
[17] A. Barrat, R. Burioni and M. Mezard, J. Phys. A29 (1996) 1311 (|cond-mat/9509142|) . 
[18] G.K.Savvidy, The system with exponentially degenerate vacuum state ( |cond-mat / 00032"2"0"| ) ■ 
[19] R.J. Glauber, J.Math.Phys. 4 (1963) 294. 

[20] G.K. Savvidy, K.G. Savvidy and PG. Savvidy, Phys.Lett. A221 (1996) 233. 



13 



